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ABSTRACT 

We analyse two-dimensional near-IR band spectra from the inner 660x315pc^ 
^ ■ of the Seyfert galaxy NGC 7582 obtained with the Gemini GNIRS IFU at a spatial 

resolution of wSOpc and spectral resolving power Rw5900. The nucleus harbors an 
^ , unresolved source well reproduced by a blackbody of temperature TwlOSOK, which 

y-^ ' we attribute to emission by circumnuclear dust located closer than 25 pc from the 

C/3 ' nucleus, with total mass of « 3 x 10"'^ M0. Surrounding the nucleus, we observe a 

ring of active star formation, apparently in the galactic plane, with radius of w 190 pc, 
an age of R:!5Myr and total mass of ionized gas of « 3 x 10^ M©. The radiation of 
\ the young stars in the ring accounts for at least 80% of the ionization observed in 

^ ' the Br7 emitting gas, the remaining being due to radiation emitted by the active 

I nucleus. The stellar kinematics was derived using the CO absorption band at 2.29 /im 

04 , and reveals: (1) a distorted rotation pattern in the radial velocity field with kinematic 

' center apparently displaced from the nuclear source by a few tens of parsecs; (2) a high 

\ velocity dispersion in the bulge of cr, = 170kms~^; (3) a partial ring of ct* = 50kms~^, 

located close to the Br7 emitting ring, but displaced by R:!50pc towards the nucleus, 
interpreted as due to stars formed from cold gas in a previous burst of star formation. 
OO , The kinematics of the ionized gas shows a similar rotation pattern to that of the stars, 

■ plus a blueshifted component with velocities ^lOOkms"^ interpreted as due to an 

outflow along the ionization cone, which was partially covered by our observations. 
The mass outflow rate in the ionized gas was estimated as Mhii ~ 0.05 M0 yr~^, 
^\ ' which is one order of magnitude larger than the accretion rate to the active galactic 

nucleus (AGN), indicating that the outflowing gas does not originate in the AGN, 
but is instead the circumnuclear gas from the host galaxy being pushed away by a 
nuclear outflow. The flux distribution and kinematics of the hot molecular gas, traced 
by the il2 X2.22fim emission line, suggests that most of this gas is in the galactic 
plane. An excess blueshift along PA« —70°, where a nuclear bar has been observed, 
can be interpreted as an inflow towards the nucleus. We thus conclude that the H2 
kinematics traces the feeding of the AGN, while the ionized gas kinematics traces its 
feedback via the outflows. An AGN-Starburst connection in the nucleus of NGC 7582 
is supported by the ratio between the mass accretion rate and the star formation rate 
in the circumnuclear region of «0.26%, which is close to the expected relation between 
the mass of the SMBH and that of the host galaxy bulge in galaxies (the Magorrian 
relation) . 

Key words: galaxies: Seyfert - infrared: galaxies - galaxies: NGC 7582 (individual) 
- galaxies: kinematics - galaxies: outflows - galaxies: starburst 



1 INTRODUCTION 

Many recent studies on the stellar population in the 
vicinity of active galactic nuclei (hereafter AGN) report 
E-mail:rogemar@ufrgs.br an excess of young to intermediate age stars in the in- 
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Figure 1. Top left: Large scale image of NGC 7582 from iHameed fc Devereuxl lll999l). Top right: HST WFPC2 image for the central 
region of NGC 7582. The green contours are for the [O III] emission form IStorchi-Bergmann et al ] l|l999l) . Bottom right: Br7 flux map 
reconstructed from the GNIRS IFU spectroscopy. Bottom left: Spectra extracted for a circular aperture of 0'.'4 radius centred at the 
positions N, A and B marked at the bottom-right panel with the emission lines and CO absorption bandhead identified. 



ner few hun dred parsecs when compared to non-active 



galax i es (e.g. Heckman et ahl 19971: Gonzalez Delgado et ahl 



2001 



200a 



1998; Schmitt & Kinney! l2000l: ICid F ernande s et a,l. 
Storchi-Bcrgmann et al.l l200ll: iKauffmanii et al.l 
Cid Fernandes et all l2005l : iDavies et al.ll2007^ ." This excess 
supports the existence of the so- called AGN-Starburst 
connection JPerrv fc Dysonl Il985l: 



198E ; iNorman fc Scovilld 11988': 'Cid Fernandes fc TerlevichI 
19951 ). which can be understood as due to the fact that 
both circumnuclear star formation and nuclear activity feed 
on gas inflowing towards the nuclei of galaxies. The gas 
infllow, if massive enough, in its way to the centre, may 
trigger circumnuclear star formation, and, reaching the 
nucleus, also the nuclear activity. Alternatively, there may 
be a delay in the triggering of the nuclear activity, which 
would be fed via mass loss from the evolving stars. The 
AGN-Starburst connection is consistent and somewhat ex- 
pected within the presently favored scenario of co-evolution 
of galaxies and their nuclear supermassive black holes 



iTerlevich fc Melnick 



(SMBH), supported by the M— tj relation dTremaiiie et al.l 
l2002l : iFerrarese fc Merritll2000l : iGebhardt et al.ll2000l ): if "the 

bulge grows via the formation of new stars (starbursts), and 
the SMBIf grows via mass accretion from its surroundings, 
one could expect a relation between the nuclear mass 
accretion rate and the star-formation rate in circumnuclear 
starbursts. 

The mechanism through which the gas makes its way 
from the circumnuclear regions of star formation - which 
in many cases have a ring-like morphology - to the AGN 
are nevertheless still elusive. In lSimSes Lopes et al.l (|2007h . 
using optical images obtained with the Hubble Space Tele- 
scope (HST) of the inner few hundred parsecs of both active 
and non-active galaxies, the authors concluded that there 
is a strong correlation between the presence of dusty struc- 
tures (spirals and filaments) and nuclear activity in early- 
type galaxies, proposing that these structures trace the ma- 
terial in its way in to feed the active nucleus. Recently, 
this was shown to be the case for at least two low-activity 
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AGN with circum nuclear rings of s tar-formation: NGC 1097 
jFathi et al.ll2006l ') and NGC 6951 (|Storchi-Bergmann et all 
In these two galaxies, the authors have observed 
streaming motions along nuclear spirals in the Hq emitting 
gas. Another simi l ar cas e is the Seyfert galaxy NGC 4051 
where iRiffel et alJ l)2008h found streaming motions along a 
nuclear spiral in molecular gas emitting the H2A2.1218^m 
line. 

Inflows towards AGN are difficult to observe, partic- 
ularly in ionized gas, due the predominance of emission 
from outflowing gas in the vicinity of the active nucleus. 
The inflow probably occurs mostly i n "cold" gas, while the 
outflows are seen in ion i zed g a s (e.g. ICrenshaw fc Kraemen 
I2OOOI : iDas et all |2005| . bood : iRiffel et al.l I2OO6I ). Studies 
aimed to investigate the kinematics and flux distributions 
of the molecular and ionized gas in the near-infrared (here- 
after near-IR) bands indeed reveal that the H2 emitting 
gas usually has a distinct kinematics and distribution 
from that of the ionized gas. The H2 kinematics seems 
to be "colder" and, together with the flux distribution, 
suggest that the molecular gas is more restricted to the 
galactic plane, while the ionized gas seems to extend to 
higher galactic latitudes and h as a "hotter" kinematics 
dStorchi- Bergman n et al.l Il999l: iRodrfguez-Ardila et al.l 
2004 Rodrig uez- Ardila. Riffel fc Pastori^ l2005al : 



Riffel et all 120061 . 120081 ). These studies suggest that 



the feeding of AGN is dominated by inflow of cold molec- 
ular gas, while the ionized gas constributes mostly to its 
feedback via the outflows. 

This work is aimed at investigating the relation between 
circumnuclear star formation and nuclear activity, as well 
as to look for the mechanisms of feeding and feedback in 
the nearby Seyfert galaxy NGC 7582, which has a Hubble 
type SBa b and h as a dista n ce of 21.6 Mpc (for z = 0.00525 
from Ide'V aucoul eurs et al.l (Il99ll ). with I'.'O corresponding 
to 105 pc at the galaxy, adopting Ho = 73kms~^ Mpc~^). 
This galaxy is an ideal candidate for t his study because it 
has both circumnuclear star formation (iRegan fc Mulchaeyl 
1999'; 'Storchi-Bcrgma nn et aT]l200ll : ISosa-Brito et al.ll200ll : 
Wold fc Galliano 2006a|)" and a nuclear outflow observed by 



Morris et all (|l985l ') using the TAURUS Fabry-Perot spec- 



trometer. The outflows, observed in the [O ill] emission line 
are oriented approximately perpendicular to the galactic 
plane, with velocities of up to — lOOkms"^, and are co- 
spatial with an ioni zation cone observed in an [O ill] A5007 
narrow band image l|Storchi-Bergmann fc Bonat to 1991). 

Besides mapping the [O ill] outflow. iMorris et al.l (119851 ) 
also reported the presence of a kpc-scal e rotating disc in 
the plane of the galaxy, observed in Ha. IWold fc Gallianol 
(2006a,) have studied the star formation in the circumnuclear 
disc using [Nell]12.8/im narrow band images, and report the 
presence of embedded star clusters claiming that there are 
no counterparts detected at optical or near-IR wavelengths. 
They calculate the mass of the SMBH as 5.5 x 10^ Mq using 
high resolution mid infrared spectroscopy. 

We use K-band spectroscopic observations obtained 
with the Integral Field Unit (IFU) of the Gemini Near- 
Infrared Spectrograph (GNIRS) to map the stellar and 
gaseous flux distributions and kinematics of the inner few 
hundred parsecs of NGC 7582. The K-band spectra allow 
the mapping of the stellar kinematics via the CO A2.3/xm ab- 
sorption feature, the ionized gas flux distribution and kine- 



matics via the Br7 emission line and the molecular gas flux 
distribution and kinematics via the H2 emission lines. This 
paper is organized as follows: in Section 2 we describe the 
observations and data reduction. In Section 3 we present the 
two-dimensional flux distributions and kinematics for both 
the emitting gas and stars. The results are discussed in Sec- 
tion 4, and in Section 5 we present our conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

The data were obt ained at the Gemi n i South Telescope with 
the GNIRS IFU l|Elias et al.l 1 19981 : lAUington-Smith et al.l 
|2006| . I2OO7I) in 2005 Nov 10 as part of program GS-2005B- 
Q-65, and comprise four individual exposures of 900 s cen- 
tred at A=2.24^m. The 111 1/mm grating with the Short 
Blue Camera (0'.'l5/pixel) was used, resulting in a nom- 
inal resolving power of R=5900. The GNIRS IFU has a 
rectangular field of view, of approximately 3'.'2 x 4'.'8, di- 
vided into 21 slices. At the detector, the slices are di- 
vided along their length into 0'.'15 square IFU elements. 
The major axis of the IFU was oriented along position 
angle PA=203°, in order to partially cover the ionization 
cone observed in the [Oiii]A5007 narrow-band image by 
IStorchi-Bergmann fc Bonattol l|l99 j ). Two sets of observa- 
tions were obtained, one centred at l'.'2 from the nucleus 
along PA=23° and another centred at 0'.'75 from the nucleus 
along PA=203°, the small offset being required to fill in the 
gaps between the slices. The observing procedure followed 
the standard Object-Sky-Sky-Object dither sequence, with 
off-source sky positions since the target is extended. Con- 
ditions during the observations were stable and clear, with 
image quality in if-band, in the range 0'.'45-0'.'55, as ob- 
tained for the full width at half maximuum (FWHM) of the 
stellar profiles. The resulting spectral resolution, measured 
from the FWHM or the emission lines of the arc calibrations 
is «2.9A. 

The data reduction was accomplished using tasks in the 
GEMINI. GNIRS IRAF package as well as generic iraf tasks. 
The reduction procedure included trimming of the images, 
flat-fielding, sky subtraction, wavelength calibration and s- 
distortion correction. The telluric bands were removed and 
the frames flux calibrated by interpolating a black body 
function to the spectrum of the telluric standard star, ob- 
served immediatelly after the galaxy. The final IFU dat- 
acube contains 840 spectra covering the spectral range from 
2.14 /im to 2.33 /xm, each spectrum corresponding to an an- 
gular coverage of 0'.' 15x0'.' 15, which translates to 16xl6pc^ 
at the galaxy. The total observed field of view 6'.'3x3'.'0 (ob- 
tained by mosaicing the two set of observations) corresponds 
to a region of projected dimensions 660pcx315pc at the 
galaxy. 



3 RESULTS 

In the top- left panel of Figure[T]we present a lar ge scale op- 
tical _R -band image of NGC 7582 from Hamccd fc Devereuxl 
l|l999h . obtained with the Cerro Tololo Inter-American 
Observatory (CTIO) 1.5 m telescope. The top-right panel 
presents a continuum image observed with the HST Wide 
Field Planetary Camera 2 (WFPC2) using the fiher F606W 
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obtained from the HST archive (program 8597 - PI: Re- 
gan, M.). We have overlaid on the HST image the con- 
tours (in green) of the narrow- band [O ill] A5007 image from 
IStorchi-Bergmann fc Bonattd (flQQl ). which maps the ion- 
ization cone. The central rectangle shows the field-of-view 
of the GNIRS IFU observations, which covers only a small 
fraction of the ionization cone observed in [O ill] . In the 
bottom-right panel we present a Br7 flux distribution im- 
age reconstructed from the GNIRS IFU datacube, where 
we have marked , in ma genta, the star clusters found by 
IWold fc Gallic { 2006al ) . These clusters delineate an elon- 
gated ring of recent star formation, with semi-major axis of 
« l'.'7, which was almost entirely covered by our observa- 
tions. 

In the bottom-left panel, we show three charateristic 
IFU spectra obtained within circular apertures of 0'.'4 ra- 
dius: the nuclear spectrum - centred on the location cor- 
responding to the peak flux in the continuum (position N 
in the Br7 intensity map), and spectra from two locations 
at 1'.'7 NW (position A), and at l'.'7 SE of the nucleus (po- 
sition B), located in the star-forming ring. The spectra of 
the latter regions present emission lines of Br7 A2.1661 /im, 
H2 A2.2235 /xm and H2 A2.2477 /im, used to map the gaseous 
kinematics and flux distributions, and the stellar absorption 
band of ^^CO (2,0) A2.2935 ^m, used to obtain the stellar 
kinematics. The nuclear spectrum is redder than those from 
the ring, with a flux in the continuum 10—15 times higher 
than in the ring, and shows a narrow Br 7 emission line, on 
the top of what seems to be a broad component of the line. 

3.1 Emission-line flux distributions 

Flux distributions in the emission lines were obtained di- 
rectly from the datacube by integrating the flux under each 
emission-line proflle and subtracting the continuum contri- 
bution using images obtained in windows adjacent to the 
emission line. 

We have also adjusted Gaussian profiles to the emission 
lines in order to measure the gas radial velocity (from the 
central wavelength of the line) and velocity dispersion (from 
the width of the line). The fitting of each emission- line profile 
was obtained using the SPLOT iraf task, and resulted also in 
an integrated flux for each line, from which we also obtained 
flux distributions. A comparison of the flux distributions ob- 
tained via the fitting of Gaussians to the flux distributions 
obtained directly from the datacube showed that, although 
both flux distributions agree with each other within the un- 
certainties, the ones obtained via Gaussian fltting are less 
noisy away from the nucleus. These flux distributions are 
shown in Figure[2]for the Br7 (top panel) and II2 A2.2235 yLxn 
(middle panel) emission lines with mean uncertainties of 6% 
and 11%, respectively. The Br7 emission is extended along 
the major axis of the galaxy (PA~ 165°) and show flux en- 
hancements in an arc shaped stru cture which approximatel y 
trace the star clusters observed bv lWold fc GaUianol (|2006al ) , 
shown again by the magenta circles. Two strong Br7 emis- 
sion peaks in the ring are located at « 1'.'7 NW and ~ 1'.'7 SE 
of the nucleus, and are encircled by the green contours in 
Fig. [2] which correspond to the limits within which the flux 
in each pixel is higher than 50% of the peak flux in the re- 
gion. We note that each of our encircled regions incl ude two 
of the clusters identified bv lWold fc GaUianol (|2006al ). While 



Bry X2.1661/tm 




Hj X2.2235Mm 




2 0-2 



Continuum 




2 0-2 
(arcsec) 

Figure 2. 2D flux maps for the Br7 (top), and H2A2.2235^m 
(middle) emission lines with mean uncertainties of 6% and 11%, 
respectively, and reconstructed image for the 2.17 /xm continuum 
emission (bottom). The black contours are iso-intensity curves, 
and the green contours in the Br7 flux map delimit regions where 
the line flux per pixel is higher than 50% of the peak flux in the 
region. The magenta filled circles indicate the position of th e em- 
bedded star clusters observed by IWold fc GaUianol ll2006al 'l. The 
black circles surrounding the nucleus in the Br7 flux map rep- 
resent the apertures used to extract the spectra shown in Fig.|8] 
and the green dashed lines show the orientation of the pseudo-slits 
used to construct one-dimensional cuts shown in Fig. 1101 



the Br7 flux distribution show intensity variations and en- 
hacements in the position corresponding to the star clusters, 
mapping the ring of star-formation, the H2 flux distribution 
(middle panel of Fig. [2} is more uniform througouth the nu- 
clear region, showing no enhancements at the star clusters. 
The black regions correspond to pixels for which the S/N 
ratio was too low to allow a reliable measurement of the 
emission line. 

The bottom panel of Figure[2] shows the 2.17 /im contin- 
uum image obtained from interpolation between two contin- 
uum windows adjacent to the Br7 profile. The iso-intensity 
contours in the continuum image are assymetric, being more 
extended to the W-SW than to the E-NE, as expected, due 
to the strong dust absorption observed in optical images, in 
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Figure 3. Fits of tlie stellar kinematics for the positions A (top) 
and B (bottom) shown in Figure[T] The observed spectra arc 
shown in black, the fitted template in red and the residuals in 
green. 

Ro od agreement with the 1.6 /im continuum image presented 
bv lWold fc GaUiand ()2006ah . 

3.2 Stellar kinematics 

In order to obtain the stellar line-of-sight velocity dis- 
tribution (LOSVD) we fitted the ^^CO(2,0) stellar ab- 
sorption bandhead at 2.2935 fj,m in the K-hand spec- 
tra using the penalized Pixel -Fitting (pPXF) method of 
ICappeUari fc EmseUemI (|2004l ). The algorithm finds the 
best fit to a galaxy spectrum by convolving template stel- 
lar spectra with a given LOSVD. This procedure outputs 
the stellar radial velocity, velocity dispersion and higher- 
order Gauss-Hermite moments for each spectrum fitted. 
The pPXF method allows the use of several stellar tem- 
plates, and also to vary the relative contribution of the 
different templates to obtain the best fit, thus minimiz- 
ing the problem of template mismatch. For this study, 
we have selected as template spectra those of the Gem- 
ini spectroscopic library of late spectral type stars ob- 
served with th e GNIRS IFU using the grating 111 1/ mm 
in the j^'-ba n d JWinge, Riffel fc Storchi-Bergmannll2007^ . In 
iRiffel et all l|2008t ). we describe some properties of this li- 
brary and discuss the problem of template mismatch using 
the CO absorption bandheads with the pPXF method. 

In Figure|3] we show the resulting fit for positions A 
(top) and B (bottom) marked in Fig.[T] For the regions where 
the signal-to- noise (hereafter S/N) ratio was lower than ~ 30 



we replaced the spectrum of each single pixel by an average 
including the 9 nearest pixels. Even adopting this procedure, 
we could obtain measurements of the stellar kinematics only 
for the central 4'.'8 x 3'.'0 region, as the S/N ratio of the 
spectra close to the borders of the IFU field becomes too 
low to obtain a reliable fit. For this same reason, we could 
not obtain stellar kinematic measurements also very close 
to the nucleus (radius<0'.'3), as the CO bandhead feature 
becomes diluted by dust and line emission. 

In Figure|4]we present 2D maps of the stellar kinemat- 
ics. Black pixels in this figure correspond to the masked 
regions where we did not obtain a reliable fit. The central 
cross marks the position of the nucleus, defined as the locus 
of the peak of the if— band continuum and the dashed line 
shows the PA=165°, which is the positi on angle of the ma- 
jor axis of the galaxy (|Morris et al.|[l985l ). The top left panel 
of Fig.|4] shows the stellar radial velocity field, from which 
we subtracted the systemic velocity of the galaxy listed in 
the HyperLedcQ database - Vs = 1605 km s"'^ (relative the 
local standard of rest), which agrees with our value for the 
nuclear region. The mean uncertainty in the radial veloc- 
ity is ~ 15kms~^. The radial velocity field shows bluehisfts 
to the SE and redshifts to the NW, suggesting rotation, al- 
though with the kinematical centre displaced from the peak 
of the continuum emission. The velocity field has an am- 
plitude of « 200 kms"^ (from -100 kms"^ to 100 kms"^) 
and seems not to be symmetric relative to the nucleus, with 
the highest blueshifts to the S-SW being reached closer to 
the nucleus than the highest redshifts to the N-NW. A one- 
dimensional cut along the galaxy major axis suggests that 
the turnover of the rotation curve to the S-SW is reached 
already at «100-150pc from the nucleus, but the turnover 
to the N-NW seems to be beyond the border of the field 
(see Fig.[lOl). 

In the top right panel of Fig.|4] we present the stellar 
velocity dispersion (a*) map, with values ranging from 40 
to 170 km s"'^, and mean uncertainties of ~ lOkms"^. The 
highest values are observed to the E and S of the nucleus. 
A partial ring of low a, values (~ 50kms~^) is observed 
at « 1" from the nucleus, surrounding it from N to SW. 
The hs and Gauss-Hermite moments maps, which mea- 
sure deviations of the LOSVD from a Gaussian distribu- 
tion, are presented at the bottom panels of Fig. |4] and have 
mean uncertainties of 0.03. The values of /13 and h4 are be- 
tween —0.18 and 0.18, which are similar to typical values 
obtained for a large number of galaxies i n similar analyses 
(jEmsellem et al.ll2004l : ICanda et aLll2006l) . 

3.3 Gas kinematics 

In the left panels of Figure[5] we present the radial veloc- 
ity fields obtained from the peak wavelengths of the Br7 
and H2 A2.2235 ^m emission lines, with mean uncertainties 
of 6kms~^ and 9kms~^, respectively. The systemic veloc- 
ity of the galaxy has been subtracted from all the emission- 
line velocity maps. The velocity field of the Br7 emitting 
gas is very similar to that of the H2 emitting gas, be- 
ing also similar to the stellar velocity field, with the S-SE 
side approaching and the N-NW side receeding. The largest 

^ http://leda.univ-lyonl.fr 
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Figure 4. 2D maps of the stellar kinematics. Top left: radial velocities with mean uncertainties of 15kms~^. Top right: velocity 
dispersions with uncertainties of lOkms"'^. Bottom: /13 and hi Gauss-Hermite moments with uncertainties of 0.03. The dashed line 
represents the PA of the major axis of the gal axy, while the dotted and so lid lines show the orientation and extent of the nuclear bar and 
nuclear spiral arms, respectively, observed by lAcosta-Pulido et al. The cross marks t he position of the nucleus and the magenta 

filled circles indicate the position of the star clusters observed bv lWbld fc Gallianol l l2006al ). The orientation is shown in the top left 
panel. 



blueshifts are observed S of the nucleus reaching values of 
~ — 170kms~^, while the largest redshifts reach smaller val- 
ues of ~ lOOkms"^ and are observed NW of the nucleus. 

The gas velocity dispersion {a) maps are shown in the 
right panels of Fig. [5] Mean uncertainties are 7kms~^ for 
Br7 and lOkms"^ for H2. The values were corrected for 
the instrumental broadening (crinst ~ ISkms"'^), and range 
from CT ft; 20 to « 90kms~^ for Br7 and from cr ft; 20 to ft; 
70kms~^ for H2. The Br7 velocity dispersion show regions 
of low values {a ~ 30kms~^) to the N of the nucleus, and 
higher values to the S-SE. The highest a values observed in 
the H2 hue are located around the nucleus, and the lowest 
values in regions away from it. 



3.4 Velocity Slices 

In order to better map the gaseous kinematics, and as 
our data has a relatively high spectral resolving power, we 
could "slice" the line profiles in a sequence of velocity chan- 
nels. Our goal is to better sample the gas kinematics over 
the whole velocity distribution, including the emissio n-line 
wings fsee lGerssen et al.ll2006l : iRiffel et al.ll2006l . l2008l ). The 
flux distributions in the velocity slices were obtained after 



subtraction of the continuum contribution determined as the 
average of the flux from both sides of the emission lines. Each 
slice corresponds to a velocity bin of « 50kms~^ (two spec- 
tral pixels) and the resulting flux distributions are shown in 
Figs. [6] and [T] In these figures, the flux levels are presented 
in logarithmic units, and the slices trace the gas from neg- 
ative (blueshifted, top left) to positive (redshifted, bottom 
right) velocities relative to the systemic velocity. For the 
Br7 emitting gas, the highest blueshifts of ft; — 360kms~^ 
are observed at SW of the nucleus, while the highest 
redshifts of ft; 260kms~^ are observed at ft; 2" NW of the 
nucleus. For the H2 emitting gas, the highest blueshifts reach 
ft; — 300kms~^ and the highest redshifts ft; 160kms~^, ob- 
served aproximately at the same locations as for the Br7 
emision. While the highest blueshifts are observed displaced 
from the major axis, in the velocity range from ft; —200 to 
200kms~^ the emission peak moves from S-SE to N-NW 
essentially parallel to the major axis, suggesting that, for 
this range of velocities, the emission is dominated by gas 
rotating in the plane of the galaxy. 
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Figure 5. Left: Radial velocity maps for the Br7 (top) and H2 A2.2235 /im (bottom) emitting gas, with mean uncertainties of 6 and 
9kms~^, respectively. Right: Velocity dispersion maps for the same emission lines with mean uncertainties of 7kms~^ for Br7, and 
lOkms"^ for H2. The dashed line represent the orientation of the major axis of the galaxy. 
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Figure 8. Nuclear spectrum of NGC 7582 for a circular aperture 
of O'.'S diameter (top spectrum), extranuclear spectrum for a ring 
with inner radius of 0'.'4 and outer radius of 0'.'6 from the nucleus 
(bottom spectrum) and the difference between the two (middle 
spectrum). The extranuclear spectrum was scaled by the ratio 
between the areas of the nuclear and extranuclear apertures. 



4 DISCUSSION 

4.1 The nuclear spectrum 

The flux distribution in the continuum shows an unresolved 
source at the nucleus (whose location we assume to coincide 
with the peak of the flux distribution) . In order to isolate the 
spectrum of the nuclear source, we have extracted its spec- 
trum within an aperture of radius 0'.'4, which corresponds to 
the radial distance where the point source profile decreases 
to 1/3 of its peak value. In order to subtract the underlying 
stellar contribution to the nuclear flux we have extracted a 
spectrum within a ring with inner and outer radii of 0'.'4 and 
C/.'6, respectively, surrounding the nuclear aperture. We show 
the apertures corresponding to the nuclear and extranuclear 
spectra as circles in the top panel of Fig.O This choice of 
apertures was dictated by the need to isolate the nuclear 
spectrum from that of the circumnuclear ring of star for- 
mation. Larger apertures for the nuclear and extranuclear 
spectra would include contribution from the star-forming 
ring. 

The resulting nuclear and extranuclear spectra are 
shown in Fig.[S] (top and bottom spectra). In the same 
Figure, we also show the difference between them (mid- 
dle spectrum) after scaling the circumnuclear spectrum to 
the flux it would have if obtained through the same aper- 
ture as that of the nuclear spectrum. We note that the 
circumnuclear spectrum does not differ much from those 
of the star-forming regions A and B, suggesting that re- 
cent star formation is occurring also inside the ring. This 
is supporte d by previous studies of the st ellar population in 
the optical dStorch i-Bergmann et al.ll200ll ) and near-infrared 
iDavies et al. 2005). 

The nuclear spectrum seems to present both a narrow 
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(arcsec) 

Figure 6. Velocity slices along the Br7 emission line profile in steps of 50kms~^. The dashed lines represent the orientation of the 
major axis of the galaxy and the velocity corresponding to the centre of each bin is shown at the top left corner of the corresponding 
panel. 



and a broad component in the Br7 emission line. The fit 
of a Gaussian to the broad component, as illustrated in the 
top spectrum of Fig.[8l gives a FWHM ^ 3900 kms"\ with 
the centre of the Une blueshifted by ^ 650kms~^ relative to 
the narrow component. The presence of a broad component 
in Br7, wi th FWHM »4000 k ms~^ has been previously re- 
ported bv lDavies et al.l |200^, in long-slit observations ob- 
tained with the Very Large Telescop e. A broad component , 
observed in Ha, was also reported bv lAretxaga et a]Tl|l999^ ■ 
although this component seems to have been produced by 
a transient event. These authors noted the appearance and 
subsequent fading of the broad Ha component as well as 
other optical emission lines and proposed three possible sce- 
narios for its origin: capture and disruption of a star by a 
SMBH, a reddening change in a dusty torus surrounding the 
nucleus, and a Type Hn supernova exploding in a compact 
starburst. 

The nuclear spectrum is redder than the extranu- 
clear one, and the subtraction of the stellar contribu- 



tion does seem to have eliminated the stellar 2.3/im 
CO absorption, clearly observed prior to the subtrac- 
tion. The resulting spectrum shows a red continuum sim- 
ilar to those of other red unresolved components ob- 
served in other AGN and attributed to emission from 
hot dust, possibly distributed in a toroidal structure 
around the SMBH, with temperature c l ose to the dust 
subli m ation t emperature (iBarvainisI Il987l: iMarco fc AUoinl 
19981. I2OOOI: iRodrfguez-Ardila. Contini fc ViegasI l2005bl : 
Rodriguez- Ardila fc Mazzalavll2006l ). Under the assumption 



that the red unresolved component is indeed due to emis- 
sion by dust, we have derived its temperature and mass, 
by fitting a blackbody function to its spectrum. The result- 
ing fit is shown as a dashed line in the middle spectrum of 
Fig. [HI corresponding to a temperature T ~ 1050 ± 140 K, 
for which the maximuum of the blackbody curve is reached 
at « 2.76 /^m. 

We have estimated the mass of the hot dust produc- 
ing the emission in the unresolved component using the for- 



-19 



-18.5 



AGN-Starburst connection in NGC7582 9 
■18 -17.5 -17 



O 

T 



^^^^^^^^^^H^ll 1 1 1 1 1 1 1 1 1 1 1 ■ 1 ■ 










S3j 






-2 2 -2 2 

(arcsec) 

Figure 7. Same as Fig.|6]for the H2 A2.2235 fim emission-line profile. 
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malism of iBarvainisI (|l987|), for dust composed by grains 
of graphite. The infrared spectral luminosity of each dust 
grain, in ergs^^ Hz^^, can be written as: 



(1) 



where a is the grain radius, Qi, is its absorption efficiency 
and Bv{Tgr) is its spectral distribution assumed to be a 
Planck function for a temperature Tgi-. 

The total number of graphite grains can be obtained 
from: 



Nhd 



(2) 



where L™ is the total luminosity of the hot dust, obtained 
by integrating the flux under the function fitted to the 
nuclear-extranuclear spectrum and adopting a distance to 
NGC7582 of d = 21.6 Mpc. In order to obtain Lf' we have 
adopted Tgr = 1050 K, as derived above, a — 0.05 / x m an d 
= 1.4 X IQ-^^iy^-^ (for u in Hz) as in IBarvainisI l|l987h . 
Finally, we can obtain the total mass of the emitting dust 
as: 



Mhd 



47r 3,^^ 

— a Nu-Dpgr 



(3) 



Adopting Pgr ~ 2.2 6 g cm ^ for the graphite density 
(|Granato fc Daneselll994l 'l. we obtain Mhd « 2.8 x 10"^ Mq. 
In Table[T]we compare this dust mass with masses of hot dust 



obtained for other AGN in previous studies. The hot dust 
mass in NGC 7582 is similar to those obtained for Mrk 766, 
NGC 1068 and NGC 3783, while it is smaUer than the values 
for Mrk 1239, NGC 7469 and Fairall 9, and higher than those 

derived for NGC 1 566 and NGC 4593. 

For Mrk 1239 jRodn'guez-Ardila fc Mazzala.vll2006l) and 
Mrk 766 (|Rodriguez-Ardila. Contini fc ViegasI bOOSbl ) the 
dust mass values were derived from long-slit spectroscopy, 
with extraction apertures of the nuclear spectrum corre- 
sponding to more than 100 pc at the galaxies, while the 
spatial resolution of our data allow us to constrain the lo- 
cation of the hot dust as being within « 25 pc from the 
nucleus. As pointed out above, the if— band excess is prob- 
ably due to the emission of hot dust from a circumnuclear 
torus which may be located even closer to the SMBH. In 
the ca se of NGC 1097, for example, IStorchi-Bergmann et al.l 
(|2005l ) have found a heavily obscured starburst closer than 
9pc from the nucleus, which could be located in the out- 
skirts of the torus. In the case of NGC 1068, interferometric 
mid-IR observations constrained the hot dust emission to 
originat e much closer, in a region of aproximately 1 pc di- 
ameter l|jaffe et al.ll2004l 'l. 



4.2 Circumnuclear star-forming regions 

NGC 7582 is known to present a kpc-scale Hq emit- 
ting disk which harbors recent star formation plus heavy 
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Table 1 . Masses of hot dust found in AGNs 



Galaxy A^jjD (1^0 ) Reference 



NGC 7582 


2.8 X 


IQ- 


3 


This work 


Mrkl239 


2.7 X 


IQ- 


2 


RodrfKucz-Ardila & Mazzalav ( 2006) 


Mrk 766 


2.1 X 


IQ- 


3 


RodriBucz-Ardila. Contini & Viegas f2005b) 


NGC 1068 


1.1 X 


IQ- 


3 


Marco & AUoin ('2000') 


NGC 7469 


5.2 X 


IQ- 


2 


Marco & AUoin (1998) 


NGC 4593 


5.0 X 


10" 


4 


Santos-Lleo et al. ( 1995") 


NGC 3783 


2.5 X 


10- 


3 


Glass (1992) 


NGC 1566 


7.0 X 


10- 


4 


Baribaud et al. (1992) 


Fairall 9 


2.0 X 


10- 


2 


Clavel, Wamsteker & Glass fl989") 



dust obscuration in the n uclear reg ion ([Morris et al.l 



Regan fc MulchaevI [l999l; fetorchi-Bergmann et al 



Ijergmann et 

Sosa-Brito et all I2001I : IWold fc Galliano! 12006^ The pres- 
ence of dust to the NE is responsible for the assymetric 
contours of the continuum flux map in the bottom panel 
of Fig. [21 The rece nt star-formation i n the dust-embedded 
clusters reported bv lWold fc GaUiand (l2006al ) (shown as ma- 
genta circles in Fig. are seen as enhancements in our Br7 
image, delineating a ring of major axis of RiS'.'S and minor 
axis of r;1'.'5. This gives an axial ratio of ^ ~ 2.3, which is 
similar to the axial ratio of the galaxy (-1 « 2.4j)) suggest- 
ing that the star-forming ring is circular, with a radius of 
«190 pc and is located in the plane of the galaxy. As pointed 
out bv lWold~fc Galliano (2_006a), these star- forming clusters 
had not been previously seen neither in the optical nor in the 
near-i nfrared. The near-IR IFU data from ISosa-Brito et al.] 
(|200ll ). although consistent with our results, does not have 
the necessary spatial resolution to resolve these clusters. 

In our data alone, clear peaks in the Br7 flux distribu- 
tion are located in the ring at l'.'7NW and l'.'7SE of the 
nucleus, and are encircled by the green contours in Fig.[2j 
which comprise the pixels with fluxes higher than 50% of the 
peak flux. We will refer to each of these two regions - whose 
physical boundaries we define as corresponding to the green 
contours - as circumnuclear star-forming regions (CNSFRs), 
and will use our data for these regions to characterize the 
star-formation in the circumnuclear ring. From the locations 
of the star clusters found bv iWold fc GaUiand (|2006al ) it can 
be seen that each of these CNSFRS correspond to two star 
clusters, thus we should bear in mind that the parameters 
we obtain for these regions correspond to two star clusters 

instead of one. 

Following iRiffel et all (|2008l ). we can estimate the mass 
of ionized hydrogen in each CNSFR as: 



1985 



2001 



) 



X 10 



Br7 



erg s~^cm" 



d 



Mpc 



(4) 



where Fsr-y is the line flux, d is the distance of the galaxy 
and we have assumed an electron temperature T — 10^ K 
and electron density = 100 cm"^ for the CNSFRs. 

We can also obtain the emission rate of ionizing photons 
for each CNSFR using (|Osterbrockll 19891 ): 



OtB Lf 



(5) 



from HyperLeda database f http:/ /leda.univ-lyonl.'frl l 



where as is the hydrogen recombination coefficient to all 
energy levels above the ground level, a'^Ha is the effec- 
tive recombination coefficient for Ha, h is the Planck's 
constant and lyHa is the frequency of the Ha line. Using 



as = 2.59 X 10"" cm'^s" 



1.17 X 10-"cm^s-i and 



Lhci/ Lbt-i ~ 102.05 for case B reco mbination assumin g an 
electronic temperature — 10* K (IOsterbrocklll989l ). we 
obtain: 



Q[H^ 



7.47 X 10 



erg s-^ 



(6) 



Finally, we can also ob tain the star formation rate 
(SFR) using (|Kennicuttll 19981 ): 



SFR 

Moyr-: 



8.2 X 10" 



L/Br'y 

ergs~i 



(7) 



under the same asumptions adoped to derive equation[S] 

As there is considerable obscuration in the nu- 
clear region of NGC 7582, the emission-line fluxes used 
in the above equations have been corrected for red- 
dening. We have used E(B — V) = 0.6, obtain ed by 
ISchmitt. Storchi-Bergmann fc Cid Fernanded (|l999l ) in a 
spectrum extracted withing a n aperture of 2x2 arcsec'^ , 
adopting the extinction law of ICardelli. Clavton fc Mathid 
(1989). 

Besides calculating the above quantities for the two 
CNSFRs, we have also calculated them for the nucleus and 
for the complete ring of star formation, as follows. In or- 
der to obtain the emission rate of ionizing photons for the 
nucleus, we integrated the Br7 emission- line flux within a 
circular aperture of 0'.'4 plus that emitted by the outflowing 
gas, under the assumption that the outflowing gas is ion- 
ized by the active nucleus. The flux in the outflowing gas 
was estimated as the sum of the fluxes in the velocity slices 
with \v\ > 200kms~^. In this manner, we obtain only a 
lower limit for the rate of ionizing photons emitted by the 
nucleus, as our field of view covers only a very small por- 
tion of the outflow, as discussed in Sec.[3l An upper limit 
can be obtained under the assumption that the total flux in 
the outfiow is the one observed times the ratio between the 
total area covered by the outflow (estimated from the [O ill] 
image) and that covered by our observations, which is ~5. 
Considering that there should be a counterpart outflow to 
the opposite side of the galaxy plane, probably hidden by 
the dust in the plane, the Br7 flux in the outflow may be 
RilO times the one measured from our data. The Br7 flux of 
the star-forming ring was obtained by integrating its value 
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Figure 9. Equivalent w idth of Br7 as pr edicted by evolutionary 
photoionization models ijDors et al.ll2008l ). 

over the whole IFU field and subtracting the contribution of 
the nuclear flux plus that of the outflowing gas. 

Tabled presents the results of the above calculations, 
together with the area used in the integrations and the Br7 
emission-line luminosities corrected for reddening, For the 
nucleus we list both the lower limit (including only the con- 
tribution of the observed outflow) and the upper limit (in- 
cluding the contribution of the total estimated outflow) . We 
also list the equivalent widths of the Br7 emission line. 

Tabled shows that the Br7 flux in inner «400 pc of 
NGC 7582 is dominated by emission from the circumnuclear 
star-forming ring, which ranges from 4 to 13 times that of 
the nucleus plus outflow. The Br7 fluxes of the two CNS- 
FRs contribute with 50 % of the total emission from the star 
formation ring. 

The emission rates of ionizing photons obtained for 
the two CNSFRs are in good agreeme nt with values ob- 
tained for CNSFRs in ot her galaxies (e.g. iGalliano &: AUoinl 
l2008l : iHagele et all 1200 7^. as well as with rates derived for 
NGC 7582 in the Mid-IR (Wold & Galliano 2006a). Such 
rates correspond to about 1000 06 stars in each CNSFR 
l|Osterbrocklll989l ). Considering that each CNSFR actually 
corresponds to two clusters, then this would correspond to 
«500 O stars per cluster. 

The star-formation rates for the CNSFRs of SFR ^ 
0.23 — 0.28 Mo/yr^^, caracterize a moderate star-forming 
re gime, and are in g ood a greement with the values derived 
bv lShi. Gu fc Pend l|200d ) for a sample of 385 CNSFRs in 
galaxies covering a range of Hubble types. 

4.2.1 Age of the CNSFRs 

We can use our data also to estima te the age of the C NSFRs, 
applying the method described bv lDors et al. I l|2008l ). which 
consists of comparing the observed values of EW(Br7) to 
those predicted by evolutionary photoionization models. As 
stressed by those authors, the real age of the CNSFRs can 



only be obtained if the contributions both from the under- 
lying bulge and from older stars in the CNSFRs are taken 
into account. 

In order to obtain the contribution of the underlying 
bulge, we used aperture photometry in the continuum image 
(Fig. El to measure both the flux of the CNSFRs continua 
and that of the bulge in the surrounding regions. We find a 
bulge contribution of «50 % of the total flux at both CNS- 
FRs, and thus the value of EW(Br7) increases from 15 A to 
30 A after this correction. 

The presence in the CNSFRs of stars formed in previ- 
ous bursts of star formation is supported by the w ork of 
ISchmitt, Storchi-Bergmann fc Cid Fernandej l|l999l ). who 
found a large spread of age in the circumnuclear star-forming 
regions. In order to ta ke this effec t into account, we apply 
the same approach as iDors et"al] (|2008h . representing the 
CNSFR stellar content by three bursts of star formation, 
with ages of 0.1, 10, and 20 Myr, and considering two values 
for the upper limit in stellar mass (Mup= 100 and 30 Mq). 
The resulting evolutionary behaviour of EW(Br7) is illus- 
trated in Figure 13 For a value of EW(Br7)«30A, we obtain 
an age of «5.0 Myr, which is similar to the v alues found 
for th e CNSFRs in NGC 6951 and NGC 1097 bv lDors et all 
(j2008t ). 

4.3 Stellar Kinematics 

As observed in Fig. |4] the stellar velocity field seems to be 
dominated by rotation. We have thus tried to fit the stellar 
velocities by a model of circular orbits in a plane subject to a 
Plummer potential, which sucessfuUy reproduced the stellar 
kinematics of the ce ntral regions of other Seyfert galaxies i n 
previous studies (e.g. lBarbosa et al.ll2006l : lRiffel et al]|2008l l. 
Nevertheless, in the case of NGC 7582 this model fails, as the 
residuals are too large everywhere. The kinematic centre of 
the rotation in Fig|3] seems to be displaced from the nucleus 
by ~ 0'.'5NW (~50pc at the galaxy), an effect which can 
also be observed in one-dimensional cuts of the velocity field 
shown in Fig. 1101 This is an interesting result which should 
be further investigated. Could this be the kinematic cen- 
tre of the galaxy? If this is the case, then the peak flux of 
the infrared continuum would not correspond to the cen- 
tre of the galaxy but instead to the location of the SMBH 
(the unresolved red source), and the ring of star formation 
would be centred on the SMBH instead of the galaxy cen- 
tre. We believe this is unlikely, as the radius of influence of 
the SMBH is much smaller than the radius of the ring, and 
thus one would expect that the ring would be centred on 
the galaxy nucleus. Our field of view is too small for us to 
reach a firm conclusion, but one possibility is that the ap- 
parent displacement of the centre of the rotation field can 
be due to distortions in a circular velo city field. Observa- 
tions in the mid-IR l|Acosta-Pulido e"t~ al. 2003) show that 
NGC 7582 presents both a nuclear bar and two nuclear spi- 
ral arms, which could be related to the observed distortions. 
The nuclear bar is oriented along PA?a 110 - 290(-70)° 
extending from ~ 1" E to ft; 1" W from the nucleus, while 
the nuclear spiral arms are oriented approximately along the 
major axis of the galaxy and extending to ft; l'.'5NW and 
SE of the nucleus. We show the orientation and extents of 
the nuclear bar and spiral as dotted and solid white lines, 
respectively in Fig.|4l We note that the SE arm seems to cor- 
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Table 2. Physical parameters obtained for the CNSFRs in NGC 7582. 



Parameter l'.'7NW I'.'TSE Nucleus'°"''+"'°^° Ring 

(nuc+llo)b ^ 



Area (arcsec^) 


1.48 


1.15 


0.5 


18.4 


logL{Br7) [ergs"-'-] 


38.53±0.03 


38.45±0.03 


37-81p^:^6j±0.04 


39.09±0.03 


log(9[//+][s-l)] 


52.40±0.03 


52.32±0.03 


51-68S^?:^^j±0.04 


52.96±0.03 


SFi?(M0yr-i) 


0.28±0.02 


0.23±0.02 




1.01±0.07 


EW(Br7) [A] 


14.7±1.4 


13.6±1.3 


O.SitO.l 


7.1±0.6 


Mhii (IQS Mq) 


8.2±0.6 


6.8±0.5 


1.5|«:^]±0.2 


29.4±2.0 



" Nucleus plus upper limit for the outflowing gas component. 
Nucleus plus lower limit for the outflowing gas component. 



respond to the location where there seems to be an "excess 
blueshift" in the stellar radial velocity field in Fig.|4l 

The stellar velocity dispersion shows the highest val- 
ues of «170kms~^ to E and S of the nucleus (Fig. [3)) 
which we interpret as representing the cr* of the galaxy 
bulge. We use this value to estimate the mass of the SMBH 
(Mbh) of NGC 7582 from the relation log(Ms///Mo) = 
Q + /31og(a*/cro), where Q = 8.13 ± 0.06, P ^ 4.02 ± 0.32 
and ao = 200 kms"^ jTremaine et al.1 |2002| ). For a, ^ 
170kms~^, we obtain Mbh ~ 7x 10^ Mq , i n good a greement 
with the value obtained by I Wold etahl (|2006bl ) through 
modelling of the kinematics of [Nell]12.8/im emitting gas. 

The (T* map (Fig|4} shows a region with low values 
(a* ~ 50kms~^) delineating a partial ring around the nu- 
cleus and which seems to be embedded in the bulge. These 
locations are clos e to tho se of the star clusters observed by 
IWold fc Gallianol (|2006al ). but somewhat displaced, being 
« 50 pc closer to the nucleus. The lower a, values reveal 
that these stars have a "colder kinematics" than that of the 
stars of the bulge. We note that these a, values are similar 
to those of the H2 gas, suggesting that these stars have re- 
cently formed in the circumnuclear ring from kinematically 
cold gas, and which have not yet "thermalized" with the 
"hotter" stars in the bulge. But in order to pre sent CO ab- 
sorpt ion, the stars should be at least 10^ yr old (jOliva et al.l 
I1995I ). what suggests that they belong to a previous burst 
of star-formation and not the one which is presently ion- 
izing the gas. This interpretation is supported by the dis- 
placement of this ring relative to that in the ionized gas, 
which suggests that the star-formation is propagating out- 
wards from the nucleus. Ring- like regions wit h low at values 
have been observed also around other AGN l|Barbosa et al.l 
[2OO6. : ,Riffel et al..,2008. ) . 



4.4 Gas Kinematics 

The radial velocity maps of both the Br7 and H2 emitting 
gas (Fig. [SJ are similar to the stellar one (Fig. |3| . In order 
to better observe the possible differences, we have extracted 
one-dimensional cuts from the radial velocities maps which 
are shown in Figure lTUI These cuts were obtained by aver- 
aging the velocities within pseudo-slits of 0'.'45 width pass- 
ing through the nucleus. The orientations of the slits were 
selected to provide a good coverage of the velocity fields, 
and are shown as dashed green lines in the bottom panel of 
Fig. [2] The top-left panel of Fig. llOl shows the velocities along 



PA= — 15(345) — 1 65°, which is the o rientation of the major 
axis of the galaxy (jMorris et al.|[l985h . while in the top-right 
panel we present the velocities along PA= 23 — 203° , in the 
bottom-left panel along PA= 60 — 240° and in the bottom- 
right panel along PA=113 - 293°. 

The highest differences between the gaseous and stel- 
lar velocities are observed at 0'.'5— 2'.'5 SW from the nu- 
cleus along PA= 203°, where the gaseous velocities reach 
values of up to — 160kms~^, while the stellar velocities 
reach at most -100kms"\ For the PAs 165° and 240° the 
gaseous velocities are still more blueshifted than the stellar 
but only by about 40kms~^ (to SE-SW of the nucleus). 
Along PA=293(-67)°, to NW-W of the nucleus, the veloc- 
ities observed for the stars are similar to those of the Br7 
emitting gas, while the H2 emitting gas is blueshifted by 
~ — 20kms^^. We note that this orientation coincides with 
that of the nuclear bar, and suggest that, as this blueshift is 
observed against the far side of the galactic plane, we could 
be observing an infiow of molecular gas towards the nucleus 
along the bar (which could be in the galaxy plane). We have 
already observed similar H2 infiows towards the nucleus of 
another active galaxy, NGC 40 51, not along a n uclear bar, 
but along a nuclear spiral arm iRiffel et al.ll2008h . 

From Fig. [To] it can be concluded that the emitting 
gas is blueshifted relative to the stars to SW of the nu- 
cleus. This region is c o-spatial with part of the [O ill] ioniza- 
tion cone observed bv lStorchi-Bergmann fc Bonattol |l993), 
shown by the green contours in Fig.[T] We interpret this ex- 
cess blueshift as due to the contribution of gas outflowing 
from the nucleus to the total gas emission. We note that we 
are observing only a small fraction of the outflowing gas be- 
cause the field-of-view of our observations is much smaller 
than the dimensio ns of the ionization cone, along which the 
gas is outflowing ( Morris et al.|[l985l ). 

Figs, m [5] and [TO] suggest that the radial velocity field 
of the gas is similar to that of the stars except for the re- 
gion along the ionization cone/nuclear outfiow and possibly 
along the nuclear bar where there may be an H2 inflow. In 
order to isolate the outflowing and inflowing components, 
we have subtracted the stellar velocity field from the Br7 
and II2 velocity fields. These residual maps are shown in 
Figure im In the residual Br7 velocity map the velocity val- 
ues to NE-N-NW of the nucleus are close to zero, while to 
the S-SW-W the gas is blueshifted by up to -lOOkms"^ 
relative t o the stars. This v elocity is of the order of those 
found bv lMorris et~aLl (|l985l ) for the [O ill] emitting gas. In 
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Figure 10. One-dimensional cuts of the gaseous (black circles 
and crosses) and stellar (open circles) radial velocities obtained 
along the PAs -15 - 165°, 23 - 203°, 60 - 240° and 113 - 293° 
for a pseudo-slit of 0'.'45 width. 



the residual H2 velocity map, the blueshift at PA ~ 293(-67) 
can be observed in approximate alignment with the nuclear 
bar, whose orientation is shown as the white solid line in the 
bottom panel of Fig. 1111 

Even higher blueshifts - of up to —350 kms~^ - are ob- 
served to the S-SW in the velocity channel maps (Fig.|S]). 
The difference between the radial velocity maps and the ve- 
locity channel maps is due to the fact that the first give the 
peak velocity of the emission line for a given spatial position, 
while the latter show how the gas with a certain velocity is 
distributed along the NLR. The peak wavelength will usu- 
ally correspond to the kinematics of the component (s) with 
the highest flux(es), while the channel maps will reveal also 
weaker components which contribute to the wings of the 
emission-line profiles. We thus attribute the high blueshifts 
seen in the velocity channel maps to low mass, high velocity 
outflows, originated close to the nucleus which compress and 
push the gas further out to the observed radial velocities of 
~ — lOOkms^^. These high velocity outflows do not appear 
in the radial velocity maps due to their low flux, appearing 
only as a blue wing to the proflle, thus not affecting the peak 
wavelenght used to derive the radial velocity. We have ob- 
served similar results in a recent st udy of the NLR of ot her 
nearby Seyfert galaxies using IFUs (|Barbosa et al.ll2009l ). 

Outflows from the nucleus of active galaxies have been 
obser ved at op tical wav elengths usin g long-slit spectroscopy 
(e.g. IStorchi-B crgman n et al.l 1 19921 : ICrenshaw fc Kraemerl 
I2OO0I; iDas et al.i i2005i. 1200^ and from radio observations 
(e.g. iMorganti et al.|[2007h and more recently by o ur group 
using IFUs l|Riffel et al.ll200^ : iBarbosa et"aLll2009l ') . 
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Figure 11. Residual maps obtained by the difference between 
the Br7 and H2 radial velocity fields and the stellar one. The 
white line in the bottom pane l show the orientation and extent 
of the nuclear bar observed bv lAcosta-Pulido et al. 1 1I2OO3I) . 



Mass outflow rates 

We have used the velocity channel maps of Fig.|S]in order 
to estimate the mass outflow rate in ionized gas. Inspect- 
ing these maps, it can be concluded that radial velocities 
of up to ±200kms~^ are observed to both sides fo the nu- 
cleus, suggesting that up to these velocities, we are mapping 
gas motions in the galaxy disk. Above 200kms~^, there are 
mostly blueshifts and almost no redshifts, what can be inter- 
preted as due to the presence of onflows. Our interpretation 
for this results is that these outflows are seen in blueshift 
to the S-SW because they are partially oriented towards 
us and are projected against the far side of the galaxy disk 
while the counterpart redshifted outflows are not seen be- 
cause they are hidden behind the near side of the dusty 
disk. Thus, assuming that the emitting gas with blueshifts 
higher than 200kms~^ is outflowing from the nucleus, we 
can estimate the mass outflow rate. In order to do this, we 
calculate the mass flux through a circular cross section with 
radius r ~ 0'.'7 at l'.'2 from the nucleus, estimated from the 
shape of the light distribution in the velocity channels cor- 
responding to —308 and — 257kms~^ in Fig.|6] The mass 
outflow rate has been calculated as Mhii ~ Ira^ v -rrr^ /, 
where is the proton mass, A^e is the electron density, 
V is the v elocity of the outflowing gas and / is the fiU - 
ing factor (|Storchi-Bergmann et al.li2007l : [Riffel et ahllioosh 
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and the factor of 2 is to take into account the likely pres- 
ence of a hidden counterpart outflow to the one we ob- 
serve. Adopting A^e = 100 cm~^ and / = 0.001 and adding 
the contribution of the two velocity bins above, we obtain 
Mhii ~ 4.8 X 10~^ Mq yr. 

This outflow rate is somewhat smaller than 
those quoted in the literature for active galax- 
ies, which are between 0.1 a.nd lOMnyr"^ (e.g. 
IVeilleuxT Cecil fc Bland-HawthornI Hool). However, the 
latter values are generally obtained for objects with higher 
levels of activity than NGC 7582. For similar levels of 
activity and in comp arable scales of hundreds of parsecs, 
iBarbosa et al.l (|2009l ) have obtained similar ionized gas 
outflow rates of « 10"^ - 10"^ M© yr"^ from 2D mapping 
of the gas kinematics in the central region of 6 nearby 
Seyfert galaxies, in good agreement with the value obtained 
here. 

We now calculate the mass accretion rate necessary to 
produce the observed nuclear luminosity of NGC 7582. The 
latter can be estimated from m — where Lboi is the 

bolometric luminosity, r] is the efficiency of conversion of the 
rest mass energy of the accreted material into radiation and 
c is the speed of light. We adopt 77 = 0.1, a typical value 
for geometrically thin and opt ically thick accretion discs 
(e.g. (Frank. King fc Rahi3l2002l) and use the relation Lboi ~ 
20 Lx (|Elvis et all 1 1994 ), w here 



{ Storchi-Bergmann et al.l 19991: Rodriguez- Ardila et al.l 



« 7.2 X 10^" erg s- 



^ ^ IS 

the hard X-ray luminosity (|Piconcelli et al.ll2007l ). The re- 
sult is rh fa 2.6 x 10"^ M© yr"\ 

Comparing the mass outflow rate with the accretion 
rate calculated above, it can be concluded that the former is 
an order of magnitude la rger than the latter, a relation sim- 
ilar to those reported bv lVeilleux. Cecil fc Bland-HawthorrJ 
l|2005l ) for m ore luminous a c tive g alaxies, as well as to those 
obtained by iBarbosa et al.l l|2009l) fo r other nearby Seyfert 

jalaxi es and to that obtained by ICrenshaw fc Kraemerl 

20071 ) for the Seyfert galaxy NGC4151. 



4-4-2 Feeding vs. feedback 

Although the are similarities between the velocity fields ob- 
served for H2 and Br7, there are also differences in the de- 
tailed kinematics for the two gas phases, as follows. (1) The 
small but systematic excess blueshift in the H2 emitting gas 
along PA= 293(— 67)°, suggests there could be an inflow of 
molecular gas along a nuclear bar, as discussed above. (2) 
The Br7 a map (Fig. [S} shows the highest values to the S- 
SE of the nucleus, approximately at the same location where 
the highest blueshifts are observed (see Fig. Ilip . This result 
supports an origin for at least part of the Br7 emission of this 
region in disturbed gas which is not confined to the galaxy 
plane but is outflowing from the nucleus along the ioniza- 
tion cone seen in [O ill]. The H2 (J map, on the other hand, 
presents overall lower a values consistent with a dominant 
origin in molecular gas orbiting in the plane of the galaxy. 
(3) The velocity channel maps shown in Figs. [6] and [7] also 
support a larger contribution from outflowing gas to the Br7 
than to the H2 emission, since the blueshifted emission in 
Br7 reaches larger velocity values and is found farther away 
from the line of nodes than in the H2 emitting gas. 

Previous studies have also found that the H2 
traces less disturbed kinematics than the ionized 
gas in the circumnuclear region of Seyfert galaxies 



20041: iRodrfguez-Ardfla. Riffel fc Pastorizal l2005al : 
iRiffel et al.ll2006l . I2OO8I ). leading to the suggestion that the 
molecular gas, orbiting in the plane, traces the feeding of 
the active nucleus via inflow of cold gas, while the ionized 
gas traces its feedback via the outflows. In I Riffel et al.l 
(2008) we could trace streaming motions in the H2 emitting 
gas towards the nucleus along a nuclear spiral arm. In the 
present observations, there is also some evidence of inffow 
along a nuclear bar. The compactness and high inclination 
of the circumnuclear ring of star-formation does not favor 
the detection of possible additional streaming motions 
between the ring and the nucleus. The II2 flux distribution 
shows that the molecular gas is not concentrated in the 
star-forming ring, but is distributed over the whole nuclear 
region, indicating the presence of similar density of molec- 
ular gas at the ring and between the nucleus and the ring. 
This molecular gas could be the material which is feeding 
and probably will continue to feed the active nucleus for 
some time in the future. 



4.4.3 AGN-Starburst connection m NGC 7582 

The total mass of hot H2 in the nuclea r region of NGC 7582 
can be estimated as l|Riffel et al.ll200§ ): 



Mh 

Mr. 



2.40 X 10' 



erg s ^ cm ^ 



/ d 



^Mpc 



(8) 



where we assume a vibrational temperature of T^n, 
2000 K, w hich implies a pop ulation fraction f^^i j^2 



3.7x10"^ l|Scovifle et al 



Asm = 2.53 X 10"'^ s"^ 



Il982l). and a transition probability 



Turner. Kirbv-Docken fc Dalgarng 



I1977I ). Integrating the flux of the H2A2.2235/im emission 
line over the wole IFU held we obtain FH2A2.2235 ~ 5.5 x 
lO-^'^ergs-^cm"^ and thus « 620 M©. 

The above value is the mass of hot H2 emitting gas, 
which is probably only the "hot skin" of the molecular gas 
available in the nuclear region to feed the active nucleus. 
The total mass of molecular gas should be much higher, as it 
should be dominated by cold gas, which has been estimated 
to be 10^ to 10^ times that of the hot gas in galaxies for 
which both hot a nd cold molecular gas have been observed 
dPale et al.ll2005l ). Thus, there should be plenty of molecular 
gas, not only to feed the AGN but also to give origin to the 
recent bursts of star-formation in the ring. 

The molecular gas accumulated in the nuclear region 
was probably the fuel which gave origin to the nuclear ring 
of star formation, as well as to the present episode of nu- 
clear activity. Although both the star formation and nuclear 
activity could have been triggered simultaneously, there is 
also the possibility that the star formation was triggered 
first, and then the mass loss from stellar winds and super- 
nova explosions in the CNSFRs have been the source of fuel 
for the SMBH. It is thus instructive to estimate the avail- 
able mass to feed the SMBH under the above assumption. 
For solar metallicity and an age of 5 Myr, the total mass 
loss for the two CNSFRs is in the range «10^ - 10^ M© 
(|Leitherer fc Heckman|[l995h . Assuming a value of 10* M©, 
and that the duration of this process is a few 10® yr, one gets 
a rate of mass loss of 10~^ M©yr~^. (Note that this value 
is ~5 times larger than the nuclear accretion rate calculated 
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above, of rh ~ 2.6 x 10"'^ yr~^, within the uncertain- 
ties of the calculations.) If this mass reaches the nucleus 
along a nuclear bar (as observed in NGC 758 2) or nuclear 
spiral arms fas observed in NGC 1097 (Fathi et all [iooil ') 
and NGC 6951 (|Storchi-Bergmann et al.ll2007l ')]. with veloc- 
ities of «20— 60kms~^, this gas would take a few 10^ yr to 
reach the nucleus. This timescale alows for the coexistence of 
the nuclear activity and star-forming regions of 5 Myr under 
the hypothesis that the nuclear activity has been triggered 
by mass accretion originated from mass loss of young stars 
formed in the circumnuclear ring. Alternatively, the trig- 
gering may have happened during a previous burst of star- 
formation in the ring. The presence of such previous star- 
burst is supported by the partial ring of low ct, observed in 
the stellar kinematics. In addition, as pointed out in Sec. 14. II 
the extranuclear spectrum does not differ significantly from 
those of the star-forming ring, suggesting that there may be 
regions of recent star formation even closer to the nucleus, 
which may be contributing to the SMBH feeding as well. 

Finally, we quantify the AGN-Starburst connection in 
NGC 7582 by comparing the mass accretion rate to the 
SMBH with the star formation rate (SFR) in the circumnu- 
clear region. As pointed out in the introduction, the M— a 
relation implies that the mass of the SMBH in galaxies grows 
in proportion to the growth of the bulge. If the bulge grows 
via star formation episodes, the mass accretion rate to the 
SMBH should be proportional to the SFR. Which should 
be the ratio between the mass accretion rate and the SFR? 
The same as the between the mass of the SMBH and the 
mass of the host galaxy bulge, which was first obtained by 
iMaeorrian et al.] ( 1998). and today known as the Magorrian 
relation (jFerrarese fc Mcrrit 2001'). which has been revised 
to a ratio of 0.1% between the mass of the SMBH and that of 
the bulge. Assuming that the total SFR can be represented 
by that of the ring, we get for NGC 7582 a ratio of 0.26% be- 
tween the mass accretion rate and the SFR, a number which 
is close to that of the Magorrian relation, considering the un- 
certainties. These uncertainties include, in particular, a pos- 
sible underestimate of the SFR in the nuclear region, as our 
observation misses the top border of the star-forming ring, 
and a possible contribution of a starburst included within 
the nuclear aperture. 



5 SUMMARY AND CONCLUSIONS 

Two-dimensional near-IR K—hand spectroscopy from the 
mner 660x315 pc^ of the Seyfert galaxy NGC 7582, obtained 
with the Gemini GNIRS IFU at a spatial resolution of « 
40 pc^ and spectral resolution of ~3 A, was used to map the 
molecular and ionized gas emission-line fiux distributions 
and kinematics as well as the stellar kinematics. The region 
covered by the observations include most of a circumnuclear 
ring of star-formation in the plane of the galaxy plus a small 
part of an outfiow extending to high galactic latitutes. 
Our main conclusions are: 

• The nucleus contains an unresolved source whose con- 
tinuum is well reproduced by a blackbody function with 
temperature T~1050K, which we atribute to emission by 
circumnuclear dust heated by the AGN, with an estimated 
mass of ~ 3 X 10""^ Mq located within 25 pc of the nucleus. 



The nuclear spectrum shows also a broad component in the 
Br7 emission line with FWHMRi3900 kms"\ 

• The Br7 flux distribution is dominated by emission 
from a circumnuclear ring of star formation, previously seen 
only in mid-IR observations, with radius of ~ 190 pc, which 
contributes with ~80-90% of the Br7 emission, the rest 
being contributed by the nuclear source. Two large star- 
forming regions with age « 5 Myr contribute with 50% of 
the Br7 emission of the ring, resulting in ~1000 06 stars per 
region and a star-forming rate of ~O.25M0 yr~^. The total 
rate of ionizing photons emitted by the ring is ~ 10^^ s~^, 
the total star- forming rate is «1 Mq yr~^ and the total mass 
of ionized gas is « 3 x 10^ M0 . 

• The H2 fiux distribution is more uniformly distributed 
over the nuclear region than the Br7 fiux distribution and 
gives a mass of hot molecular gas of « 620 Mq . From previ- 
ous studies, it can be concluded that the total mass of H2 
gas (which is dominated by cold H2) can be lO^-lO'^ times 
larger. 

• The stellar velocity field shows a distorted rotation pat- 
tern whose centre appears to be displaced up to 50 pc from 
the nucleus (identified as the peak of the continuum emis- 
sion). The distortion could be associated to a nuclear bar 
and a nuclear spiral previously seen in mid-infrared images. 
The velocity dispersion of the bulge is ct* ~ 170 kms"'^, im- 
plying in a mass of Mbh ~ 7 x 10^ Mq for the SMBH, in 
good agreement with previous determinations. Immersed in 
this bulge, there is a partial ring of lower a, (~ 50kms^^), 
which is close to the star-forming ring but displaced ~ 50 pc 
inwards. As the velocity dispersion is similar to that of the 
molecular gas, this ring can be interpreted as due to stars 
recently formed from cold gas, although in a previous burst 
(age ^ 10^ yr), still keeping the gas kinematics as they did 
not have time yet to "thermalize" with the stars of the bulge. 

• The radial velocity field of the ionized gas is simi- 
lar to that of the stellar component to the N-NW, while 
to the S-SW there is an additional blueshifted component 
{v fa -lOOkms"^), which we attribute to outflows along 
the ionization cone (previously observed in a narrow-band 
[O III] image) which is partially covered by our observations. 
Velocity channel maps along the Br7 emission-line profile 
show even higher blueshifts, reaching v « — 300kms~^ 
in the outflow, while the maximum redshift observed is 
V « 200 km s^^. The velocity dispersion is enhanced in the 
region of the outflow. 

• The mass outfiow rate in the ionized gas is estimated 
to be A/hii ~ 0.05 Mq yr~^, which is an order of magnitude 
larger than the accretion rate to feed the AGN, a ratio which 
is similar to those found in other AGN. This large ratio 
indicates that the outfiowing gas does not originate in the 
AGN, but is instead the circumnuclear gas from the host 
galaxy being pushed away by a nuclear outfiow. 

• The kinematics of the hot molecular gas, traced by the 
H2 emission, shows smaller blueshifts along the outfiow as 
well as lower velocity dispersions, which suggests that most 
of the molecular gas is in the galactic plane. An excess 
blueshift along PA« —70°, where a nuclear bar has been 
observed, can be interpreted as an infiow towards the nu- 
cleus. We thus conclude that the H2 kinematics traces the 
feeding of the AGN, while the ionized gas kinematics traces 
its feedback via the outfiows. 

• The estimated mass loss rate from the evolving stars 
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in the circumnuclear ring is « 10^^ Moyr^^, which is «5 
times the nuclear accretion rate. As the estimated time for 
this mass to reach the nucleus is a few 10^ yr, this material 
could be the fuel which just triggered the nuclear activity. 
Nevertheless, as there is also a previous burst of star forma- 
tion almost co-spatial with the ring, the fuel may have been 
available previously (e.g. w lO'^yr). In addition, previous 
data suggest that there may be also recent star formation 
inside the ring, closer to the nucleus. 

• Wo conclude that the AGN-Starburst connection in 
NGC 7582 may have occurred in two ways: (1) there is a 
large molecular gas reservoir in the nuclear region which 
gave origin to both the star-formation in the circumnuclear 
ring and nuclear activity; (2) the molecular gas has given 
origin first to circumnuclear star formation and the mass 
loss from the evolving stars have then triggered the nuclear 
activity. 

• The ratio between the mass accretion rate and the SFR 
in the circumnuclear region, of 0.26% is close to the value 
implied by the Magorrian relation between the mass of the 
SMBH and the mass of the bulges of the host galaxies. This 
result shows that a growth of the SMBH proportional to 
the growth of the bulge can proceed via circumnuclear star 
formation around AGN. 
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